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ABSTRACT 

We use the Millennium simulation to probe the correlation between cluster veloci- 
ties and their shapes and the consequences for measurements of the kinetic Sunyaev- 
Zeldovich (kSZ) effect. Halos are generally prolate ellipsoids with orientations that are 
correlated with those of nearby halos. We measure the mean streaming velocities of 
halos along the lines that separate them, demonstrating that the peculiar velocities 
and the long axes of halos tend to be somewhat aligned, especially for the most mas- 
sive halos. Since the kSZ effect is proportional to the line-of-sight velocity and the 
optical depth of the cluster, the alignment results in a strong enhancement of the kSZ 
signature in clusters moving along the line of sight. This effect has not been taken into 
account in many analyses of kSZ signatures. 
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1 INTRODUCTION 

Observations of large-scale flows of matter in the universe 
provide strong constraints on structure formation theories 
and cosmology. In particular, observations of peculiar ve- 
locities of cluster s of galaxies could constrain cos mological 
parameters fe.g.. iBhattacharva fc Kosowskvll2008l . and ref- 
erences therein) or poss ibly challenge the ACDM paradigm 
iKashlinskv et aLllioOsI ). 

The new generation of Cosmic Microwave Background 
(CMB) expe riments such as the Atacama Cosmology Tele- 
scope (ACT, lKosowskvll2003l ) and the South Pole Telescope 
(SPT, iRuhl et al ll2004( ) wiU identify thousands of galaxy 
clusters out to redshifts beyo nd z ~ 1 via the therma l 
Sunyaev-Zeldovich (tSZ) effect (|Sunvaev fc Zeldovichlll972h . 
The tSZ effect in clusters arises from the inverse Compton 
scattering of CMB photons by free electrons in the intra- 
cluster medium (ICM) and involves a change of frequency of 
CMB photons, with radio photons typically being shifted to 
frequencies above 217 GHz . Another smaller effect, the ki- 
netic Sunyaev-Zeldovich effect (kSZ), which is due to the pe- 
culiar velocity of electrons in the ICM produces a frequency 
independent distortion to the CMB spectrum. Attempts to 
measure the kSZ effect have mostly produced upper limits 
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on the peculiar veloci ty of clusters l|Holzapfel et al.l Il997l : 
iBenson fc et al.] |2003| ) but the current CMB experiments 
should detect the kSZ effect with high significance (e.g., 
iHernandez-Monteagudo et al.ll200^ ). 

Prospects for extracting peculiar velocities of clus- 
ters detected via the kSZ, and using these for cosmol- 
ogy, have been discussed b y many authors: a nalyt i- 
cal work includes that b y Rephaeli fc Lahav lll99ll). 
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The potential to cross-correlate kSZ signatures with other 
foreground t r acers is discussed in iDore et al.i (2004, ) and 
iDeDeo et al.l (|2005l). We also note that , despite cl aims to 
the contrary (|Croft fc Efstathioul [l99i : iPeell |2006| ). linear 
theory should describe cluster velocity correlations fairly 
well if the mean streami ng motions of m assive halos towards 
each other are included (|Sheth fc Zehavi.2009. ). 

Almost all of the work on the extraction of kSZ signa- 
tures of clusters assumes spherical clusters but cluster ha- 
los are not perfectly spherical. In this paper we work out 
the correlation between the orientation of ellipsoidal galaxy 
clusters and their peculiar velocity in order to investigate 
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the enhancement of the kSZ signal due to this correlation. 
The kSZ signal is proportional to th e line-of-sight velocity 
and t he optical depth of the cluster (jSunvaev fc Zeldoviclil 
1 19721 ). On average an alignment between the orientation and 
velocity will result in a higher optical depth in clusters mov- 
ing along the line of sight. 

Section [2] provides a brief description of the data uti- 
lized in our calculations. Section [3] presents the correlation 
between the velocity and shape of the clusters and section |3] 
describes the effect this correlation has on measurements 
of the kSZ signal. Section [5] provides some discussion and 
conclusions drawn from our results. 



Table 1. The velocity parallel to the orientation of the cluster 
divided by the velocity perpendicular to the cluster for the various 
mass ranges. 





Mass range( h ^Mq) 




1012 _ 1013 10" -10" 10"- < 1015 








1.051 1.072 1.100 



2 SIMULATION AND HALO CATALOG 

Our h alo catalog was extra cted from the Millennium Simu- 
lation l|Springel et al.|[2005h which adopted roughly concor- 
dance values for the parameters of a flat A cold dark matter 
(ACDM) cosmological model, r2dm = 0.205 and fib = 0.045 
for the current densities in CDM and baryons, h = 0.73 
for the present dimensionless value of the Hubble constant, 
(Tg = 0.9 for the rms linear mass fluctuation in a sphere of 
radius 8 ^~^Mpc extrapolated to 2: = 0, and n = 1 for the 
slope of the primordial fluctuation spectrum. 

The halos are found using a two-step procedure. The 
first step entails identifying all collapsed halos with at 
least 20 particles using a friends-of-friends (FoF) group- 
finder with linking parameter b = 0.2 times the mean 
particle separation. T he substructure algorithm SUBFIND 
l|Springel et al.ll200ll ) is then used to subdivide each FoF- 
halo into a set of self-bound sub-halos. In this study we 
exclusively focus on sub-halos which we simply refer to as 
halos. 

We present results for halos within mass bins ranging 
from 10^2 to 10" /i-^Mq, lO" to lO" h'^MQ and lO" to 
lO^'^ ft-^Mo. There were 401334, 43139 and 2701 halos in 
each respective mass bin. 



3 SHAPE- VELOCITY ALIGNMENT 

Two general properties of clusters are well known. Firstly 
it has been shown that in general clusters tend to point 
towards each other (e.g., iFaltenbacher et"al] |2002| ) and 
secondly that they t end to stream towards each other 
(jSheth fc Zehavill2009l ). From these properties it can be an- 
ticipated that there should therefore be some correlation 
between the velocity of a cluster and its orientation. 

To investigate such a correlation we extract these quan- 
tities from the halo catalog. The velocities of the halos are 
calculated by taking the peculiar bulk velocity of all parti- 
cles belonging to the halo. As is standard, the orientations 
are derived by evaluating the eigenvalues of the inertial ten- 
sor, or more precisely, the second moment of the mass tensor 
of the halo. The maximum eigenvalue obtained corresponds 
to the major axis of the halo and the minimum to the minor 
axis. 

To confirm the correlation between halo orientations 
and neighbouring halos we compute 9, which is the angle 
between the major axis of a halo and the the connecting 
line to another. We then bin the angles derived for each 



pair according to r, which is the separation between the two 
halos. Finally we determine the mean cosine for logarithmic 
separation bins and calculate the mean angle for various sep- 
arations of halos. The results are shown in Figure[T]for halos 
within the indicated mass range. A mean cosine > 0.5 in- 
dicates alignment between connecting line and orientation. 
For separations of < 50 h~^Mpc halos tend to point towards 
each other and the amplitude of this alignment is clearly en- 
hanced for more massive halos. 

To investigate the extent to which clusters move along 
the line of separation we determine the average velocity of 
halos along the connecting line as a function of separation. 
We determine the velocities of each halo along the connect- 
ing line by first calculating the direction of the line joining 
the two halos and then calculating the velocity component 
along this connecting line. These velocities are then binned 
logarithmically according to separation to obtain the mean 
streaming motion of the halos. The results are shown in 
Figure [21 Halo pairs within the considered mass and separa- 
tion ranges tend to move towards each other with increasing 
velocities as the distance between them de creases. This re- 
sult co mpares well with results published bv lSheth fc Zehavil 
(|2009l ). 

Figure [3] displays the probability distribution of the co- 
sine of <j) where (j) denotes the angle between the orientation 
of the halo and its bulk velocity. This was calculated by 
binning each of the halos according to its cosine of (j). For 
all mass ranges we find an excess at large values of cos{(j}) 
indicating a preference for small angles between orienta- 
tion and velocity. We see a correspondingly small number 
of halos with large angles between their orientations and ve- 
locities. A significant deviation from isotropy, displayed by 
the dotted horizontal line, for all mass ranges is apparent. 
Although this deviation increases with halo mass, in gen- 
eral the orientations of halos above lO^'^ h~^MQ tend to be 
aligned with their velocities. We note that observations of 
the SZ effect are mostly only possible in clusters with masses 
M > lo^* h~^MQ where the effect is most pronounced. For 
this reason, in the following section where we investigate the 
kSZ effect, we do so only for the most massive halos. 

The extent to which the velocities are aligned with the 
orientation has been quantified in Table[T]where the average 
velocity parallel to the major axis for each mass bin is calcu- 
lated and divided by the average velocity perpendicular to 
the major axis. One can see that for larger halos the veloc- 
ity parallel to the major axis is on average 10% higher than 
the velocity perpendicular. This trend is smaller for smaller 
mass halos as their elongation is less dramatic. 
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Figure 1. Mean cosine of angles, 9, between halo orientations 
and the connecting line to neighbouring halos, for halos within 
the indicated mass range, as a function of separation. The 
dotted line illustrates the results for an isotropic distribution. 



Figure 2. Mean streaming velocities of halos along the line 
of separation within the indicated mass ranges. 
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Figure 3. Probability distribution for the cosine of the angle, 
0, between the orientation and the velocity of halos within the 
indicated mass ranges. An isotropic distribution would result in 
the dotted horizontal line. Poissonian errors are shown. 

4 IMPACT ON KSZ MEASUREMENTS 

The kSZ effect is proportional to the optical depth of a halo, 
T, and the velocity of the halo along the line of sight, . The 
optical depth is defined by 

T = O-T [ ne{l)dl ^ [ Alpdm{l) 

J nip S2o J 

where or is the Thomson scattering cross section, is the 
electron number density along the line of sight and rUp is the 
m ass of the proton. In defining this we assume, as outlined 
in iDiaferio et al.1 (|2000| ). that the electron number density 
is related to the dark matter mass density along the line of 
sight, pdm.il), by the following equation 

pdm(l) Slfc 

n-e ~ ^7^- 

TTlp iio 

where ili, is the baryonic density and Qq is the mean matter 
density in units of the critical density. We have thus assumed 
that the gas traces the dark matter halos and relates the 



properties of those halos to the properties of the gas. It 
is important to note that the fraction of gas at the cores 
of clusters may be lower than the cosm ic value by 30% 
l|Kravtsov et al.ll2005l : lEttori et al.ll2006l ). As a resuh of this, 
the amplitude of the kSZ signal as calculated below could 
also be reduced by this percentage. However, the relative 
ratios between the dispersion parallel and perpendicular to 
the line of sight and between the real and random samples 
as used in our investigation will not be affected. 

Peculiar velocities along the line of sight, , will pro- 
duce kSZ temperature fluctuations. This kinematic effect is 
given by 




where (l) is the bulk velocity along the line of sight at I 
and c is the speed of light. 

We calculate ksz for halos in our simulation by placing 
a box of size 5.5 /i~^Mpc centered on each halo and divid- 
ing it into cells of various le ngths. The spatial re solution of 
the simulation is 5 /i~^kpc ijSpringel et al.l[2005h while the 
best angular resolution possible with current telescopes ('^ 1 
arcminute) corresponds to about 40 /i~^kpc in the nearest 
clusters. Therefore we have investigated a range of pixel sizes 
from 50 /i^^kpc upwards. The box size and the number of 
cells per dimension were chosen to ensure that the center of 
halo coincides with the center of the most central cell. The 
integral is calculated by summing the contributions from 
discrete cells along one axis which corresponds to the line 
of sight, I. Hence we obtain 2D maps of the kSZ signal for 
each halo. Figure |4] shows an example of the kSZ signal of a 
cluster-sized halo from two viewpoints. The maximum value 
of the kSZ signal in the map, which in all cases comes from 
the central pixel, is taken as the amplitude of the kSZ signal 
for that halo and, to relate this to work on the power spec- 
trum of kSZ fluctuations, we calculate cTksz, the root mean 
square fluctuation in this amplitude. This ap proach is simi- 
lar in nature to the peak model presented bv lYoshida et all 
(|200ll ). Using their arguments we are thus confident that 
possible contributions resulting from secondary structures 
such as filaments may be neglected. 
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minor axis [Mpc/li] major axis [Mpc/li] 

Figure 4. Two projected kSZ maps of the same halo of 1.2 X 10^* /i^^Mq. Its pecuUar velocity is 640kms~^, the velocity components 
along the major, intermediate and minor axes are 597 kms^^, 57kms~^ and — 224 kms~^, respectively. The left panel shows a kSZ 
map if the halo is rotated such that the major axis is parallel to the line of sight. The right panel displays the map for the same halo 
with the minor axis oriented along the line of sight. For illustrative purposes we here use much higher resolution than that used for the 
determination of cr^sz shown in Figure [5] 



The kSZ signature of a cluster elongated along the line 
of sight is enhanced relative to a spherical cluster of the same 
mass because it has a larger optical depth. To separate this 
effect from the additional enhancement due to the velocity 
correlation, we calculated the kSZ velocity dispersion firstly 
for the correct halo peculiar velocities, and secondly for the 
case where the directions of the velocities are randomized. 
The transformation from real to random bulk velocities has 
to be performed on a particle by particle basis. First, the 
bulk velocity is subtracted from each particle's velocity and 
after a random rotation it is added to each of the particles 
again. These new velocities are used in the integration along 
the line of sight. The kSZ signal is present ed in units of 
hiK/K] for comparison with previous work (jYoshida et al.l 
I2OOII ). 

The results for two different pixel sizes are illustrated 
in Figure [5] where the coloured lines show the kSZ velocity 
dispersion for the true velocities of the halos, as a function 
of the angle between halo orientation and the line of sight. 
The peak signal is enhanced by up to 60% for halos where 
the orientation is along the line of sight compared to those 
oriented perpendicular to the line of sight. The black line 
in Figure [5] represents the kSZ signal when the cluster ve- 
locities are randomized. We observe that randomizing the 
velocities diminishes the peak signal by approximately 20%. 
The observed behavior indicates that the dependence of the 
kSZ signal dispersion on the orientation of the clusters rel- 
ative to the line of sight is due to both the change in opti- 
cal depth (if the cluster is resolved) and the velocity-shape 
alignment. The lower panel of Figure [5] displays the residual 
AfJksz/o-ksz- 

To see the effect of resolution we have computed the kSZ 
signal for additional cell/pixel sizes, including 50, 1000 and 
2000 h~^kpc on a side. For L = 50 /i~^kpc, we find higher 
peaks in the kSZ signal and an enhanced dispersion due to 
the larger optical depth for the central pixels. However, the 
dependence on orientation relative to the line of sight and 
the relative difi'erence to the cluster sample with random- 
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Figure 5. Upper panel: the rms fluctuations on a 200 /i-lkpc and 
500 /i~^kpc scale in the kSZ signal as a function of the cosine of 
the angle between halo orientation and the line of sight. Results 
are plotted for only the largest mass range {M > IQ-*^* /i~-'^Mq). 
The coloured lines show the (Tusz obtained using the velocities 
measured in the simulation. The black lines show the cr^sz ob- 
tained when the bulk velocities are randomized. Lower panel: the 
difference between the measured kSZ signal and the randomized 
signal (normalized by the randomized signal) , illustrating the ef- 
fect of the velocity-shape alignment 

ized velocities does not change. If the cell size used in the 
integration along the line of sight is large enough, we expect 
to include most of the dark matter halo and the velocities 
associated with it, thereby excluding effects due to advanta- 
geous orientation. We investigated the effect for pixel sizes of 
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L = 1 and 2 h~^Mpc and confirmed that although the am- 
phtude of the kSZ signal was reduced and the dependence 
on orientation was slightly weakened the relative difference 
to the randomized velocity sample, Aa/a, did not change 
significantly. For L — 2 h'^Mpc pixel sizes we observed only 
20% difference between the kSZ dispersions for clusters ori- 
ented parallel and perpendicular to the line of sight. For 
the corresponding randomized velocity sample the kSZ dis- 
persion is largely independent of orientation. Since at that 
resolution the entire cluster is covered by one pixel and the 
optical depth is independent of cluster's orientation, we can 
be confident that the observed dependence of the kSZ dis- 
persion on the orientation relative to the line of sight is due 
only to the velocity-shape alignment. 



5 CONCLUSION 

Clusters of galaxies have been shown to both point towards 
each other and have a tendency to stream towards each 
other. Motivated by this we have investigated the correlation 
between the orientation of a cluster and its mean peculiar 
velocity. We have found in general that this correlation is 
significant. For massive clusters we have observed that the 
velocity parallel to the major axis is on average 10% higher 
than velocity perpendicular to the orientation of the cluster 
(Table [T]). 

We have then explored how the orientation of clusters 
effects the kSZ signal and found that there is a 60% en- 
hancement of the kSZ signal dispersion for massive clus- 
ters which are orientated along our line of sight compared 
to those orientated perpendicular to the line of sight. This 
value is based on a pixel size of L = 200 ft-^kpc. The dif- 
ference for unresolved clusters, L = 2 h~^Mpc, is 20%. A 
smaller, but significant, effect for less massive clusters can 
also be seen. 

The enhancement of the kSZ signal is a result of the 
alignment of both the orientation and the velocity with the 
line of sight. In order to disentangle these effects and quan- 
tify the effect of the velocity alone we introduced a random 
sample. In this sample the orientations of the clusters re- 
mained fixed and their bulk velocities were randomized. We 
observe a 20% difference in the kSZ dispersion for the ran- 
dom sample which is caused by the dependence of the optical 
depth of the cluster relative to the line of sight alone, indi- 
cating that the remaining 20% of the enhancement is due to 
the velocity-shape alignment. 

We have observed a correlation between the orientation 
and the velocities of clusters on scales up to 100 /i~^Mpc. 
Such large-scale correlations, particularly in small-scale sur- 
veys, may result in a large bias in the kSZ signal depending 
on the orientation of large scale structures relative to the 
line of sight. In addition, since we have shown up to a 60% 
increase in the kSZ signal in a line-of-sight cluster, the se- 
lection of kSZ sources may be substantially biased towards 
clusters aligned along the line of sight. 
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